Dissipation-managed soliton in a quasi-one-dimensional Bose-Einstein condensate 
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We use the time-dependent mean-field Gross-Pitaevskii equation to study the formation of a 
dynamically-stabilized dissipation-managed bright soliton in a quasi-one-dimensional Bose-Einstein 
condensate (BEC). Because of three-body recombination of bosonic atoms to molecules, atoms are 
lost (dissipated) from a BEC. Such dissipation leads to the decay of a BEC soliton. We demonstrate 
by a perturbation procedure that an alimentation of atoms from an external source to the BEC 
may compensate for the dissipation loss and lead to a dynamically-stabilized soliton. The result of 
the analytical perturbation method is in excellent agreement with mean-field numerics. It seems 
possible to obtain such a dynamically-stabilized BEC soliton without dissipation in laboratory. 

PACS numbers: 03.75.Lm, 05.45.Yv 
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Solitons are solutions of wave equation where localiza- 
tion is obtained due to a nonlinear interaction. Solitons 
have been noted in optics p], high-energy physics and 
water waves ID, and also in Bose-Einstein condensates 
(BEC) 0, S H II 0. The bright solitons of BEC rep- 
resent local maxima [H, 0, E l m , whereas dark soli- 
tons represent local minima 0- Lj- lj] In addition to the 
observation of an isolated bright soliton in an expulsive 
potential [(| , a number of bright solitons forming a soli- 
ton train was observed by Strecker et al. |5| , where they 
suddenly turned a repulsive BEC of 7 Li atoms attractive 
by manipulating an external magnetic field near a Fesh- 
bach resonance 01 • Consequently, the BEC collapsed, 
exploded and generated a soliton train which was stud- 
ied in detail. Also, a bright vector soliton in a repulsive 
BEC su ppo rted by interspecies attraction has been stud- 
ied [HEi- 

A soliton or solitary wave by definition propagates 
over large time intervals without visible modification of 
shape which makes it of special interest. However, a soli- 
ton of BEC, or a BEC in general, suffer loss of atoms 
due to three- body recombination leading to formation of 
molecules [I3,[l3,Il3- This means that a BEC soliton will 
decay and eventually disappear as it propagates. It would 
be of interest if an artificial situation could be created in 
laboratory with a supply of atoms so as to compensate 
for the three-body recombination loss of a BEC soliton 
to generate a dynamically-stabilized soliton. To the best 
of our knowledge we demonstrate for the first time, us- 
ing the mean-field Gross-Pitaevskii (GP) equation [Tr^| . 
that such a dynamically stabilized soliton could indeed 
be prepared in a radially trapped and axially free BEC. 
As a strict soliton appears only in one dimension, we shall 
be concerned in this paper with a quasi-one-dimensional 
BEC soliton in a cigar-shaped trap in a axially symmetric 
configuration. 

To demonstrate the presence of a dynamically- 



stabilized dissipation-managed BEC soliton, we employ 
both time-dependent and time-independent analytic per- 
turbation techniques and a complete numerical solution 
of the GP equation. The numerical result is found to 
be in excellent agreement with that obtained from the 
perturbation techniques. 

Bright solitons are really eigenfunctions of the one- 
dimensional nonlinear Schrodinger (NLS) equation. 
However, the experimental realization of bright solitons 
in trapped attractive cigar-shaped BECs has been possi- 
ble under strong transverse binding which, in the case 
of weak or no axial binding, simulates the ideal one- 
dimensional situation for the formation of bright solitons. 
The dimensionless NLS equation in the attractive or self- 
focusing case pj 



sustains the following bright soliton pj: 



(0.1) 



<fi(y, t) = a sech[a(y — vt)] x 
x exp[ivy — i(v 2 — a 2 )t/2 + to], 



(0.2) 
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with three parameters. 

The parameter a represents the amplitude as well as 
pulse width, v represents velocity, the parameter o is a 
phase constant. The bright soliton profile is easily rec- 
ognized for v — as UjAy,t)\ = a sech(ay). There have 
been experimental 0, |y] and theoretical 0, S studies 
of the formation of bright solitons in a BEC. In view of 
this, here we study for the first time the possibility of a 
dynamically-stabilized dissipation-managed bright BEC 
soliton in a quasi-one-dimensional configuration. 

In recent times there have been routine experimental 
studies on the formation of BEC in the presence of a peri- 
odic axial optical-lattice potential 0, [Lsl Il9l | formed by 
a polarized standing-wave laser beam. This leads to a dif- 
ferent condition of trapping from the harmonic trap and 
generates a BEC of distinct modulation. Hence we also 
consider in this paper the modulations of a dissipation- 
managed bright BEC soliton in the presence of an optical- 
lattice potential. 
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The time-dependent Bose-Einstein condensate wave 
function ^(r, r) at position r and time r may be de- 
scribed by the following mean-field nonlinear GP equa- 
tion 111 lis 
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+ V(r) 
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-K z n 2 



*(r,r) =0, 



(0.3) 



with normalization J dr|W(r, r)| 2 = TV. Here m is the 
mass and N the number of bosonic atoms in the con- 
densate, n = \^/(r,T)\ 2 is the boson probability den- 
sity, g = Anh 2 a/m the strength of inter- atomic at- 
traction, with —a the atomic scattering length. The 
trap potential with axial symmetry may be written as 
V(r) = ^moj 2 (p 2 + v 2 z 2 ) where oj and vui are the angu- 
lar frequencies in the radial (p) and axial (z) directions 
with v the anisotropy parameter. The term K3 denotes 
the three-body recombination loss-rate coefficient and T 
denotes the constant alimentation of atoms from an ex- 
ternal source. 

For the study of bright solitons we shall reduce Eq. 
(|0.3|) to a minimal one-dimensional form under the ac- 
tion of stronger radial trapping. The one-dimensional 
form is appropriate for studying bright solitons in the 
so-called cigar-shaped quasi-one-dimensional geometry 
where v « 1. For radially-bound and axially-free soli- 
tons we eventually set v = 0. 

For v — 0, Eq. 1)0. 3|) can be reduced to an effective 
one-dimensional form by considering solutions of the type 
*(r,r) = $(z,t)V' ( ' H/ !J ) where 



\^{p)\ 2 ee ^fexp 



muip 



(0.4) 



The expression (|0.4f) corresponds to the ground state 
wave function in the radial variable in the absence of 
nonlinear interactions and satisfies 

_ ^V^ (0) + ^W 0) = M/> (0) , (0-5) 



with normalization 
2tt 



\^\p)\ 2 pdp=l. 



Now the dynamics is carried by $(z,r) and the radial 
dependence is frozen in the ground state tp^\p). 

Averaging over the radial mode if)^(p), i.e., multiply- 
ing Eq. (|0.3(l by -i/A°)*(p) and integrating over p, we ob- 
tain the following one-dimensional dynamical equation 

m. 



d h 2 d 2 muj 



+ i- 



r h m 2 ui 2 



2 37r 2 r 



1*1- 



$(z,r) = 0, 



(0.6) 



In Eq. (|0.6|) the normalization is given by 
\<f>(z,T)\ 2 dz = N 



and we have set the anisotropy parameter v = Q to re- 
move the axial trap and thus to generate an axially-free 
quasi-one-dimensional soliton. 

For calculational purpose it is convenient to reduce 
the set (|0.6f) to dimensionless form by introducing con- 
venient dimensionless variables. For this purpose we 
consider the dimensionless variables t ~ tlo, y = z/l, 
4> = y/(2aj<f>, 7 = T/(?uj), and f = K 3 /(24:n 2 a 2 l 4 uj) with 



/ = \Jfij (ujm), so that 
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<P(y,t)=ie(4 ) )<J ) (y,t), (0.7) 



where 



(0.8) 



and the normalization is given by 



\<f>(y,t)\ 2 dy = 



2aN 



In addition to studying an axially-free dissipation- 
managed soliton, we also consider a dissipation-managed 
soliton in a periodic optical-lattice potential V(y) = 



Vosin (2-Ky/X) in the axial direction jl7l Il8t UJfl, so that 
the NLS equation of interest becomes 



.d_ 



2dy 2 



<j>{y,t)=it{<l>)<i>{y,tUQS) 



where Vq is the strength and A is the wavelength of 
the laser used to generate the optical-lattice potential. 
Equation (|0.9fl will generate a dissipation-managed soli- 
ton with an axial periodic modulation. 

In Eq. (|0.7fl e(<f>) is a small perturbation. In the 
absence of a perturbation (e = 0) the soliton profile is 
known. The fate of the soliton when e ^ can be ob- 
tained by a time-dependent perturbation technique Q 
which assumes that the functional form of the soliton re- 
mains unchanged in the presence of a small perturbation. 
However, the soliton parameters change with propaga- 
tion. The most general form of the perturbed soliton is 
taken as p| 



<j>(y,t) = a(t)sech[a(t)(y - q(t))} exp[ia(t) - iyS(t)]0.10) 

In the absence of a perturbation, a and S are constants 
but q(t) and a(t) arc determined by the following ordi- 
nary differential equations pj: 



dq 

~dt 



-5. 



da 
~dt 



-(a 2 -S 2 ). 



(o.n) 
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The use of the perturbation technique leads to the fol- 
lowing ordinary differential equations for the soliton pa- 
rameters a and S @|: 



da 
~dt 



dS 

— = — Im 

dt 



Re / e(</>)F{y)dy, 



(0.12) 



i<f>) tsnh[a{y-q)W{y)dy. (0.13) 



In this work we shall be concerned only with the profile 
of the soliton with zero velocity: q(t) = 0. This can be 
obtained by solvi ng Eq . ^T^i. W ith the e of Eq. (|0~%|) 
and the </> of Eq. (|Q~Ttl . Eq. (|0~T2"jl becomes 



da 
7F 



7 



16 ^ 5 

7a fa 5 

' 15 ^ 



m & dy (0.14) 
(0.15) 



We are interested in finding the constant amplitude a 
of a dynamically-stabilized dissipation-managed soliton 
for large t from the solution of Eq. (|0. 151) . It is assumed 
that the perturbation is switched on at t = 0. The initial 
soliton at t = is taken to satisfy a(t = 0) = ao- With 
this condition Eq. l|0.15|l can be integrated to yield 



In 



16 1 4 

15 { ° - 1 



-1/4 



= 7*, 



the solution of which is 



15dg7e 



lOa^e 4 ^ - (16^ - 15 7 ) ' 



(0.16) 



(0.17) 



which for finite £ and 7 leads for large time t to the 
amplitude 



15 7 



"I 1/4 



16£ 



(0.18) 



The most interesting feature of this result is that the 
lowest-order perturbation-theory prediction for the am- 
plitude of the dissipation-managed soliton is independent 
of the initial choice for the amplitude ao at t = 0. The 
dissipation-managed soliton is robust and depends only 
on the ratio 7 /£ of the dissipation parameters 7 and £ and 
is independent of the initial choice ao in the lowest-order 
perturbation theory. 

The above result that a robust time-independent soli- 
ton of amplitude H0.18|) can be formed at large times in 
the presence of a weak dissipation suggests that the same 
could be derivable from a time-independent perturbation 
analysis. For a "stationary state" <p(y,t) = 4>(y) exp(if/,t), 
from Eq. (|0.7|) one could write the following time- 
independent equation 



id* 

2V 



= ie(M(i/), (0-19) 




FIG. 1: (Color online) The dissipation-managed dynamically- 
stabilized soliton |</>(t/)| °f the NLS equation 10.71 vs y at 
different times for £ = 0.01, 7 = 0.0002 and a = 1/2. The 
initial soliton at t — is taken as <f>(y) = 0.5sech(a;/2). The 
displayed perturbed solution labelled 'perturb' is in excellent 
agreement with the numerical result. 



where fi is a chemical potential. Employing the usual 
time-independent perturbation theory, for Eq. (|0.19l) to 
have the soliton l|0.10fl with q(t) — as a solution, one 
should have 



(M<f>)\4>) 



e(0)|0(y)| 2 rfy = 0. 



(0.20) 



With e{<t>) given by Eq. I|0.8fl and <j>{y) — a sech(ay), 
the solution of Eq. (|fT20"|) is given by Eq. flHHfr . This 
demonstrates the equivalence between the results of time- 
dependent and time- independent perturbation theories. 
However, it should be recalled that perturbative result 
(|0.18|) is valid only for small dissipation e((f>) or for small 
values of 7 and £. 

We solve the NLS equation l|0.7|) for dynamically 
stabilized dissipation-managed bright solitons numeri- 
cally using a time-iteration method based on the Crank- 
Nicholson discretization scheme elaborated in Ref. p2l 
|23| . We discretize the mean- field- hydrodynamic equation 
using time step 0.01 and space step 0.1. 

We performed the time evolution of Eq. I|0.7|) start- 
ing with the solution <fi(y) = 0.5sech(y/2) exp(£i/8) with 
ao = 1/2 for e(<f>) = at t = 0. The dissipation e{4>) 
is introduced for t > and we look for the dissipation- 
managed soliton for large time t. If a converged solution 
is obtained, it corresponds to the dissipation-managed 
soliton. We also repeated the calculation with a different 
initial choice <f>(y) = sech(y) exp(ii/2) with ao = 1 to ver- 
ify if the dissipation-managed soliton is independent of 
the initial choice as predicted by the perturbation treat- 
ment above. 

In the first numerical simulation we take £ = 0.01, 7 = 
0.0002 and ao = 1/2. The t = solution is taken as 
4>{y) = 0.5sech(y/2). Upon a time evolution of the NLS 
equation i|0.7fl the dissipation-managed soliton emerges 
at large times. In Fig. 1 we plot the numerically- 
calculated dissipation-managed soliton at times t = 
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FIG. 2: (Color online) Same as in Fig. 1 but with £ = 0.002 
and 7 = 0.00025. 



2000,4000,6000, and 8000 as well as the result of the 
lowest-order perturbation theory. At large times the 
numerically-calculated dissipation-managed soliton ex- 
hibits excellent convergence properties and agrees very 
well with the result of the lowest-order perturbation the- 
ory. We also repeated this calculation with a different 
<zq(= 1). In agreement with the result of perturbation 
theory, the dissipation-managed soliton for ao = 1 is in- 
distinguishable from that for ao = 1/2. Although the 
dissipation-managed soliton appears like a stable solu- 
tion of the NLS equation l|0.7|) , it cannot be called a sta- 
tionary solution of a time-independent Schrodinger-type 
equation with real eigenvalue. It is rather a dynamically- 
stabilized solution of the dissipative NLS equation (|0.7|) 
where there is a delicate cancellation between the two 
imaginary dissipative terms £ and 7 to maintain a stable 
profile of the soliton. 

To study how the dissipation-managed soliton changes 
with a change of the dissipation parameters £ and 7 we 
repeat the calculation of Fig. 1 first with smaller values 
of these parameters: £ = 0.002 and 7 = 0.00025. A very 
similar scenario emerges for the dynamically-stabilized 
soliton which we exhibit in Fig. 2. However, in this case 
the dynamically-stabilized soliton has a larger amplitude 
a(= 0.58) governed by Eq. I|0.18fl than in Fig. 1 where 
a(= 0.37). In these numerical studies we verified that 
for small dissipation the dissipation-managed soliton is 
independent of ao and is solely determined by the ratio 
7/£. However, for large dissipation parameters this is 
not quite so and quite expectedly the numerical solution 
of the soliton does not agree with the result of lowest- 
order perturbation theory. Also, with the increase of 
dissipation parameters it becomes increasingly difficult 
numerically to find a stabilized soliton. With further 
increase in dissipation parameters we could not obtain 
a stabilized soliton. Numerically, no stabilized soliton 
could be found for £ = 0.1. 

Finally, we consider a dissipation-managed bright soli- 
ton formed on a periodic optical-lattice potential con- 
trolled by Eq. (|0.9fl . In our simulation we take Vq = 5 



FIG. 3: (Color online) The dissipation-managed soliton of 
Eq. J£U in the presence of an optical-lattice potential with 
Vo = 5 and A = 2tt. The other parameters are the same as in 
Fig. 1. 



and A = 2tt. To solve Eq. I|0.9|) with this optical-lattice 
potential we take the initial soliton at t = to be the 
one calculated in Fig. 1 above for V(y) = and consider 
the time evolution of Eq. (|0.9|) . During this time evolu- 
tion the optical-lattice potential is slowly introduced in 
an interval of time t of about 2000 units, so that a sta- 
ble dissipation-managed soliton of Eq. 1|0.9I) is obtained 
at large time. We study the stability of this dissipation- 
managed soliton by continuing the time evolution. Again 
a dynamically-stabilized dissipation-managed soliton is 
obtained which is independent of the input guess for ao. 

Now let us see to what values of the three-body re- 
combination rate K3 the parameters £ = 0.01 and 0.002 
considered in this paper correspond to for typical exper- 
imental values of a = 10 nm, I = 1 /im and to = 2ir x 50 
Hz. Recalling that K 3 — 247r 2 £a 2 Z 4 w we find that for 
£ = 0.01, K 3 = 7.5 x 10~ 26 cm 6 /s and for £ = 0.002, 
K a = 1.5 x 10~ 26 cm 6 /s. Typical estimates of K 3 for 
different atoms 0, 0, 0] are comparable to or less 
than these values. Hence the present demonstration that 
dissipation-managed soliton can be obtained for typi- 
cal values of £ less than 0.01 is compatible with experi- 
mental values of K 3 . This leads to the conclusion that 
dissipation-managed solitons can be generated in labora- 
tory for realistic values of the parameters. 

To summarize, we have demonstrated the possibility 
of the generation of a dynamically-stabilized dissipation- 
managed robust soliton in a quasi-one-dimensional BEC 
under a typical experimental situation. This can be used 
in laboratory to generate robust BEC solitons where the 
loss of atoms due to three-body recombination could be 
compensated by alimentation of atoms from an external 
source. 
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